In [1]:

    

import numpy as np
import os, sys
path = os.path.abspath('../library/')
if path not in sys.path:
    sys.path.append(path)
from convert_to_recarray import create_gdr2mock_mag_limited_survey
import ebf, subprocess


    

In [2]:

    

### change files
'''
default start = 0
default end = 10
0-6 == thin disk
7 = thick disk
8 = halo
9 = bulge
what can be changed:
/norm_thin_disk
/norm_thick_disk
/norm_halo
/norm_bulge
'''
for thick_norm in [1.0,1.2,1.4]:
    print(thick_norm)
    print('#######################################')
    modelFile = '/home/rybizki/Programme/GalaxiaData/Model/population_parameters_BGM_update.ebf'
    ebf.update_ind(modelFile, '/popidstart', 7, ind=0)
    ebf.update_ind(modelFile, '/popidend', 8, ind=0)
    ebf.update_ind(modelFile, '/norm_thick_disk', thick_norm, ind=0)
    ### make new bhtree
    args = ['galaxia', '-s', 'warp']
    p = subprocess.Popen(args, stdout=subprocess.PIPE, stderr=subprocess.PIPE)	
    print("Galaxia spawns new BHtree")
    (output, err) = p.communicate()
    if True:
        print(output)

    ### calculate mock catalogue
    name = 'thick_disk_norm_%.1f' %(thick_norm)
    create_gdr2mock_mag_limited_survey(nside = 512, outputDir = '../output/%s_100pc' %(name),
                                       use_previous = False, delete_ebf = True,
                                      fSample = 1, make_likelihood_asessment=False, r_max = 0.1,
                                      verbose = False)
    create_gdr2mock_mag_limited_survey(nside = 512, outputDir = '../output/%s_0.001' %(name),
                                       use_previous = False, delete_ebf = True,
                                      fSample = 0.001, make_likelihood_asessment=False)


    

    




1.0
#######################################
Galaxia spawns new BHtree
b'Galaxia-v0.81\nCODEDATAPATH=/home/rybizki/Programme/GalaxiaData/\nReading tabulated values from file- /home/rybizki/Programme/GalaxiaData/Model/vcirc.dat\nfeh=-0.17 sig_feh=0.27\nThick disc sigma_v=[67, 51, 42, 0.33, 0.33]\n    [feh, dfeh, age, dage]=[-0.48, 0.3, 1.1e+10, 1e+09]\nThickDisk:\nMaking Tree\nAge=1.1e+10 lAge=0\n Total Nodes=318553 leafs=278734\nCompleted %  <1..10..20..30..40..50..60..70..80..90..100..>\nMass=1.16564e+10 Leafs=278734 Total Nodes=318553 1.16685e+10 20372.9\nWriting tree to file:               /home/rybizki/Programme/GalaxiaData/BHTree-2.3/bhtree_with_wf/bhtree_7_E0.ebf....Done\nTime Tree generation/reading =      12.7056     \nTotal Time=                         12.7337     \n'






    




---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-2-0d1f3ccb6263> in <module>()
     24 
     25     ### calculate mock catalogue
---> 26     name = 'thick_disk_norm_%,1f' %(thick_norm)
     27     create_gdr2mock_mag_limited_survey(nside = 512, outputDir = '../output/%s_100pc' %(name),
     28                                        use_previous = False, delete_ebf = True,

ValueError: unsupported format character ',' (0x2c) at index 17

In [ ]: